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We show that the linear-stability analysis of the birth of Faraday waves on the surface of a fluid 
is simplified considerably when the fluid container is driven by a triangle waveform rather than by 
a sine wave. The calculation is simple enough to use in an undergraduate course on fluid dynamics 
or nonlinear dynamics, ft is also an attractive starting point for a nonlinear analysis. 
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I. INTRODUCTION 

Parametrically excited surface waves, or Faraday 
waves, are generated when a container of fluid is ver- 
tically vibrated. Above an acceleration threshold, waves 
appear on the fluid surface. The waves vibrate at half the 
driving frequency, a feature associated with the paramet- 
ric resonance that is responsible for pumping energy into 
the surface-wave modes. 

In most "Faraday experiments," the container is driven 
sinusoidally. That is, its vertical position may be de- 
scribed by x(t) = (xpp/2) coswt, and its acceleration by 
a(t) = —uj 2 (x P p/2)cosujt. In this paper, we explore the 
consequences of driving the container with a triangle- 
wave forcing. Although the difference between the two 
waveforms matters little physically, we shall show that it 
simplifies the analysis remarkably. 

There are two broad motivations for our study. First, 
although the theory of hydrodynamic stability occu- 
pies an important conceptual place in current courses 
in hydrodynamics, standard examples such as thermal 
convection (Rayleigh-Benard Convection, or RBC) and 
shear flow between rotating concentric cylinders (Taylor- 
Vortex Flow, or TVF) are too complicated to present 
in standard undergraduate courses. Textbooks for such 
courses (e.g., usually limit themselves to heuristic 
arguments. By contrast, our analysis is simple enough 
to present in a single class to third- or fourth-year stu- 
dents. (The calculation's difficulty is equivalent to solv- 
ing the Kronig-Pcnncy model for energy bands in a one- 
dimensional solid. This example is commonly presented 
in undergraduate solid-state-physics courses.) 

A second motivation comes from the increasingly 
prominent role that the Faraday experiment occupies as 
an example of a non-equilibrium pattern-forming system. 
Most of the work on pattern formation done to date has 
focused on convection (RBC) and shear flow (TVF) 



However, a recent series of papers has established a num- 
ber of advantages for the Faraday experiment J^-|i| . Very 
large cells (> I0 2 wavelengths wide) may be created, and 
the basic ("fast") time scale can be very short (to ;$ 1CP 2 
sec.) The short time scale is particularly important be- 
cause the dynamics of pattern-forming instabilities is of- 
ten slow. Current work looks at features that occur on 
time scales of 10 4 io to 10 5 io- Convection, by contrast, 
usually has to « 1 sec, or slower. 

As we shall show below, our variant of the Faraday 
experiment clarifies our understanding of the instability 
itself. The analysis separates into two pieces: one is the 
derivation of the surface-wave dispersion relation (Sec- 
tion III); the other is the description of how parametric 
pumping injects energy into the fluid (Section II). By 
choosing triangle-wave forcing, we greatly simplify the 
latter piece. We hope that this will ultimately allow the 
nonlinear analysis (particularly of high-viscosity fluids) 
to be pursued further than it has been hitherto. For 
this reason, we make no restriction on the fluid viscos- 
ity; the only approximation made is in linearizing the 
Navicr-Stokes equations. 

The body of this paper is organized as follows: In Sec- 
tion |l| we study the related problem of the Mathieu 
equation with delta-function forcing. Solving this prob- 
lem turns out to be mathematically equivalent to solving 
for the acceleration thresholds of Faraday waves. In Sec- 
tion [II , we apply this analysis to the Faraday experiment 
with triangle- wave forcing. In Section IV, we consider 



a simple extension, driving by an asymmetric triangle 
wave, which models the recent "two-frequency" experi- 
ments that have attracted much attention || . Section |v| 
contains a brief conclusion. 



*e-mail: bech@chaos.phys.sfu.ca 
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II. THE MATHIEU EQUATION WITH 
DELTA-FUNCTION FORCING 



An infinitely extended continuous system such as a 
fluid can be thought of as containing a continuum num- 
ber of modes. To linear order, the modes are independent 
and the system's behavior can be determined from the 
study of each mode separately. Nonlinear effects imply 
mode coupling. In this section, we examine the effects of 
parametric forcing on a single mode. In Section III, 



we 

apply the results derived here to the continuum of modes 
present in the Faraday experiment. 

There are two ways to pump energy into an oscilla- 
tor. The first is by direct forcing, where the prototypical 
equation is 



x + 2-fi + luqX = f(t), 



(1) 



where f(t) is the driving force and is often chosen to be 
/ cos ut. For small damping 7, a directly forced oscilla- 
tor has a resonant response (of amplitude ss f/ujQ^) at 

LO ~ Ll>q. 

The second way to pump energy into an oscillator is 
by parametric forcing, whose prototypical equation is the 
(damped) Mathieu equation: 



(a) 



(b) 



(c) 



FIG. 1. Driving function f(t). (a) Periodic train of delta 
functions of equal signs, (b) Periodic train of delta functions 
of alternating signs, derived from a symmetric triangle wave, 
(c) Periodic train of delta functions of alternating sign, de- 
rived from an asymmetric triangle wave. 



A. Mathieu equation driven by a periodic series of 
delta functions 



+ 2 7 ± + a;g[l + /(t)]a: = 0. 



(2) 



When f(t) is sinuosoidal, there is an infinite series of 
resonances centered on lo — 2lo$, loq, (2/3)o;o, lo/2, 
(2/n)LOo, for integer n. The strongest of these res- 
onances is the first (n = 1), and a response at half the 
driving frequency is good evidence that a system is being 
driven parametrically. Below, we shall find it convenient 
to rephrase the resonance conditions in terms of periods: 
At = (n/2)T, where lo = 2n/At and lo q = 2n/T. 

Although f(t) may be sinusoidal, the solutions x(t) 
are not. Let f(t) = /cos2oj i and assume x(t) 
,7V cosioot- The term f(t)x(t) in Eq. (^) is proportional to 
cosLOotcos2LOot and can be rewritten as (1/2) (cos tu^t + 
cos3ojo^)j implying that x(t) has a 3wot component. 
Adding this to x(t) generates a bu^t term, and so on. 
Eq. (0) must in general be solved numerically. For special 
choices of f(t), however, Eq. (||) may be solved analyt- 
ically. The simplest case is f(t) equal to a sum of delta 
functions M|. A slightly more complicated case, f(t) 
equal to a square wave, was discussed in these pages some 
years ago |n]]. In Section III, we show that analyzing a 
fluid container driven by a triangle wave is equivalent 
to studying the Mathieu equation driven by a series of 
equal-amplitude delta functions, with alternating signs. 
(See Fig. [ljb.) Here, we review the analysis of M] for the 
particular case of delta-function forcing with equal and 
then with alternating signs. 



In units where time, damping rate, and forcing are all 
scaled by loq, we wish to solve 



x + 2-yx + [1 + e S(t - nAt)] x = 0. 



(3) 



The forcing function is shown in Fig. |T|a. In between 
delta-function "kicks," the system is just a free, damped 
oscillator. Let x n (t) be the solution valid between the 
kicks at nAt and (n + l)At. One finds 

Xn {t) = Ane^-^+c.c. (4) 



with ip = — 7±iw and lu' q = \J\ — j 2 . The kick imposes 
two conditions linking x n and x„+i: 

• x is continuous. 

• the velocity x jumps discontinuously. This condi- 



over a small 



tion can be derived by integrating Eq. 
interval of time centered on t n+ \. 
The conditions may be expressed as: 

x n +i(t n +i) = x n (t„+i) (5a) 

x n +i(t n +i) - in(Wi) + £x n (t n+ i) = 0. (5b) 

Imposing the continuity and velocity-jump conditions al- 
lows one to relate A„ +1 to A n : 



A r 
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n+l 



C 

S + aC C 



S 

- aS 



4 r 

v 



(6) 



where A r n and A % n are the real and imaginary parts of A 
respectively, C = cosui' Q At, S - 
is the scaled forcing amplitude 



sina> Ai, and a — s/lj'q 
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Eqs. (||) and (||) constitute a solution to Eq. (||). 
Given initial conditions Aq, we can find x(t) by first iter- 
ating the map e _7A *M to find the proper A n . Here, M 
is the 2x2 matrix in Eq. (^|). Then use Eq. (f§ to find 
x(t). 

Eq. (^|) can also be used to find the acceleration thresh- 
old for resonant response. Parametric resonance differs 
from ordinary resonance in that it is a true instability. 
If e is below a threshold s c , x(t) — > at long times. For 
e > £ c , x(t) grows without bounds. (Adding a nonlinear 
term in x to Eq. (||) will limit the oscillator's ampli- 
tude. In practice, nonlinearities are always present in 
the physical situation being modeled.) Our goal now is 
to find £ c {At). 

The key observation is that the long-time amplitude 
of the motion x(t) is controlled by the eigenvalues of the 
matrix M: if their magnitude is greater than e 7At , the 
motion's amplitude will grow each period. If the magni- 
tude of both eigenvalues is less than e 7At , damping will 
reduce the motion each cycle and the system will tend 
towards x(t) = 0. Loosely speaking, M accounts for the 
energy pumped into the system during each period At, 
while e~ 7At represents the energy lost during At. 

The matrix M has unit determinant. If we write M in 
the general form [ (j jj \ ^ ne eigenvalues are given by 



Ai, 2 = x ± Vx 2 - 1 where x = (A + D)/2 = (l/2)TrM. 
Instability arises when both eigenvalues are real and the 
magnitude of the largest exceeds e 7At . Simple algebraic 
manipulations then give for the threshold condition 



x = ± cosh 7 At 



(7) 



In the present case, the matrix elements are A = coso> At 
and D = cosu>' At— sinw At. Inserting these into the 

threshold condition Eq. (^), we find 



£ c (At) = 2u' 



cos u>' At ± cosh 7 At 



sin Wo At 



(8) 



This threshold condition is plotted in Fig. |^a for differ- 
ent choices of the damping 7. The curves generated by 
choosing the plus sign are centered on At = (n/2)T, n — 
1, 3, 5, ... and are known as subharmonic "tongues," while 
the curves generated by choosing the minus sign are cen- 
tered on At = mT, m = 1,2,3,... and are known as 
harmonic tongues. These curves are qualitatively simi- 
lar to those of the Mathieu Eq. with sinusoidal forcing 
PPI. 
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FIG. 2. Resonance tongues for the Mathieu eq. driven by a 
periodic sequence of delta functions, for 7 = 0.0, 0.1, and 0.2. 

(a) Delta functions of equal sign, (b) Delta functions of alter- 
nating signs, (c) Delta functions derived from an asymmetric 
triangle wave, for 7 = 0.2 and different values of r = At2/At. 
For r = 0.5, the tongues are the same as the dashed lines in 

(b) . Note how the harmonic tongues are more prominent as 
the asymmetry is increased (r — * 0). 



B. Mathieu equation driven by a sequence of delta 
functions of alternating sign 

In our version of the Faraday experiment, the con- 
tainer's vertical position x{t) is a triangle wave. Taking 
two time derivatives of x(t), we see that the acceleration 
a(t) is a periodic sequence of delta functions of equal 
magnitude but alternating sign. For this reason, we gen- 



eralize the analysis of section [I A to a driving term of 
the form 



£ I P(* - nAt) -S[t-( 

{ n 



n + 1/2) At}} } x(t) (9) 



Let x n (t) be the solution valid over (nAt) < t < 
[(n+ l/2)At] and let x n+1 / 2 (t) be the solution valid over 
[(n + 1/2) At] < t < [(n + l)At]. Except for the kicks 
at (nAt) and [(n + 1/2) At], the oscillator is free and we 
have 



c n (t) = A n e^-^ + 



c.c. 



+ 1 ' 2 ' + c.c. 



(10a) 
(10b) 



Applying the continuity and jump conditions at (nAt) 
and [(n + 1/2) At], we find 
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n+1 
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n+1 



_ p -~(At/2 



ac c — as 



A 

A< 



n+l/2 _ e -lAt/2 
-1/2 / 



— s 
ac c + as 




small 7/wo, which gives At « tt/w [1 - (l/2)(7/w ) 2 ] 
( ll£ ahd w = 2w [l + (l/2)(7/w ) 2 ]. To lowest order in 7/wo, 
the resonance conditions are just the zero- viscosity ones: 
it = tt/luq and oj = 2wo- At this order, the overall forc- 
lg threshold is given by e* — ttj. 



(Hb 



where c = cos ui' At/ 2 and s = sin uj' At/ 2. Com- 
bining these, we can write A n+ i — e~^ At M2MiA n = 
e -jAt MAn _ Since detM = (- det Mi)(det M 2 ) = 1, we 

have the same threshold condition, x = ± cosh 7 At, lead- 
ing to 



e c (At) = 2w 01 



/ cos ui' At ± cosh 7 At 

1 — COSLu' At 



(12) 



The resonance tongues are shown in Fig. ||b. There are 
only subharmonic tongues, a feature that turns out to be 
special to the use of our particular driving function Eq. 

To understand more intuitively why there is resonance 
for the subharmonic case but not for the harmonic case, 
consider the two time series shown in Fig. ||. Fig ||a 
shows a time series for (At = T/2, u> = 2luq), where 
T = 2-k/ujq is the natural period of the oscillator. Fig. 
||b shows the case (At = T, ui — uq). The delta- function 
impulses are indicated by up and down triangles on each 
curve while the force vector for each kick is indicated by 
an arrow above or below the triangle. Because the delta 
functions are multiplied by the coordinate x, the sign of 
the velocity increment depends on the position of the par- 
ticle relative to the origin. In the subharmonic case, the 
impulse always increases the speed of the motion. In the 
harmonic case, the up triangles decrease the speed, while 
the down triangles add to it. Because the overall motion 
is damped, it is clear that the net effect of the kicks is to 
decrease the energy of the system. One can understand 
this better by considering the matrices M\ and Mi in the 
limit of zero damping (harmonic case). The coswQAt/2 
and sincjQAt/2 terms are then —1 and 0, respectively, 
and the product M%Mi is just the identity matrix, inde- 
pendent of the value of the forcing a. Thus, no energy is 
added to the system even as it is forced arbitrarily hard. 
(One kick adds energy while the next one removes it.) 
When damping is added, as in Fig. ||b, the system is al- 
ways stable about x = 0. The harmonic case differs from 
the harmonic case of sinusoidal driving, where there is 
still a resonance because the forcing is distributed over 
the whole period of the oscillator. 

Note that when 7 > 0, higher order tongues have 
higher thresholds than the lu — 2uio tongue. The bot- 
tom of this first tongue defines an overall threshold that 
represents the global-minimum amplitude s* required to 
excite oscillations. We can find e* (7,0/0) by differenti- 
ating Eq. (|l|) for e(At). This yields cotWoAt/2 = 
7/u/q tanh7At/2 for At*(7,Wo)- This transcendental 
equation can be solved numerically or perturbatively for 
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FIG. 3. Time series x(t) for the Mathieu eq. driven by a 
periodic sequence of delta functions of alternating signs. In 
both time series, the damping 7 = 0.1 and the forcing e = 0.7. 
The triangles show where the forcing is applied. They point 
up for positive delta functions and down for negative delta 
functions. The arrows above the triangles indicate the di- 
rection of the force vector for each kick. Because the delta 
function forcing term multiplies the coordinate x, the sign of 
the force vector opposes that of the delta function for x < 0. 
Time is in units of the undamped-oscillator period (2n). (a) 
Subharmonic case (b) Harmonic case. 



III. FARADAY WAVES WITH TRIANGLE- WAVE 
FORCING 

In this section, we shall apply the above results to the 
linear-stability analysis of the flat surface of a viscous 
fluid held in a laterally infinite container of fluid of infi- 
nite depth. (The finite-depth case is straightforward, but 
the algebra is longer |^q| .) Already in 1883, Rayleigh had 
suggested that the instability studied fifty years earlier by 
Faraday might be due to the phenomenon of parametric 
resonance, and he wrote down and explored the proper- 
ties of Eq. independently of Mathieu |L4| . Seventy 
years later, Benjamin and Ursell 112] showed that the nor- 
mal modes of the linearized, inviscid Navier Stokes equa- 
tions each obeyed a Mathieu equation. In the Benjamin- 
Ursell calculation, the effects of damping must be added 
perturbatively (using 7 = 2vk 2 /u> < 1 as the bulk dis- 
sipation coefficient, where v = r\/ p is the kinematic vis- 
cosity and k the magnitude of the wavenumber of the 
disturbance). Recently (1994) Kumar and Tuckerman 
performed a complete linear analysis that is valid for all 
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viscosities |T]|. Since the analysis for arbitrary viscos- 
ity is only slightly more complicated than that for small 
viscosity, we present the more general result. 

We begin by deriving the equations of motion, follow- 
ing Kumar and Tuckerman. The Navier-Stokes equations 
of motion for the fluid layer are 

[dt + U -VIC/ = —VP + vV 2 U - g(t)z (13a) 
P 

V ■ U = 0. (13b) 

Here, U(x,y,z,t) is the fluid velocity field in the lab 
frame of reference and P(x,y, z,t) is the pressure field. 
There is a simple solution, corresponding to a fluid at rest 
in its own frame of reference: U = and P(t) = —pg(t)z. 

Linearizing about the at-rest solution, we find for the 
perturbing fields u and p 



d t u = — -Vp - 
P 

-- 0. 



vV 2 u 



V 



(14a) 
(14b) 



(In the new frame of reference, gravity does not ap- 
pear in the bulk equations of motion. It will resur- 
face in the boundary conditions.) Applying the operator 
(z ■ V x Vx) and using the incompressibility condition 
V • u = 0, we eliminate the pressure and obtain 



(di 



vV 2 )V 2 u = 0, 



(15) 



where u = u z (x, y, z, t) is the z-component of the velocity 
field. 

At the bottom of the cell, no-slip boundary conditions 
imply that u — 0. Using the incompressibility condition, 
we rewrite this as 



d z u = 0. 



(16) 



At the air-fluid interface, there are three conditions. 
The first is a kinematic one that relates the fluid velocity 
to the position of the interface: 



[d t + u ■ V]C = u\z=c, 



(17) 



where ((x,y,t) is the location of the interface. To lin- 
earize, we drop the u-V term and evaluate the right-hand 
side at z = 0, giving <9 t £ — u\ z= q. 



The second condition asserts that there are no tangen- 
tial stresses at the interface. Since the interface normal 
is just z to first order, one has that n xz — -K yz = at 
z = 0, where the stress tensor 7r is defined as 



V (diU 3 



d 3 Ui) 



-p + pg(t)z] S i: j + r\ (diUj 
(18) 



with i and j ranging over x, y, and z, and where r\ is the 
ordinary fluid viscosity. Taking the horizontal divergence 
of Eq. (18), we have 



V 2 H u + d z {V H ■ u H ) = (V| - d zz )u = 0, 



(19) 



where ur is the horizontal velocity field at z = and 
the incompressibility condition was used to write the last 
identity. 

The third condition relates the normal stress to the 
Laplace pressure caused by surface tension acting on a 
curved interface: 



7« 



-7VfrC 



(20) 



where k is the mean curvature of the interface and the 
second term comes from writing the curvature to first- 
order in £. We have, at z = 0, 



p(0 = 2 V (d z u)+pg(t)( + 1 V 2 H (. 



(21) 



Taking the horizontal divergence of the Navier-Stokes 
equation, one can obtain a second relation involving the 
pressure field and use this to rewrite Eq. (^T|) in terms 
of u and C alone: 



(d t - vV 2 )(d z u) = 2isV%(d z u) + g(t)V 2 H ( - ^V%(- 

P 



(22) 



Because the horizontal eigenfunctions of a laterally in- 
finite container are e ±jfc r , we can replace V|j by —A; 2 , 
and V 2 by (d zz — k 2 ). We then have the following lin- 
earized equations of motion for u: 



[d t - v{d zz - k 2 )] (d zz - k 2 )u = 

u = d z u = 
(0 ZZ + k 2 )u\ z=a =0 

[dt - v(d zz - k 2 ) + 2vk 2 } (d z u)\ z=0 = 

dtC = u\ z =q 



z = — oo 
z = 



g[l - f(t)\ - Ik 3 



(23a) 
(23b) 
(23c) 

(23d) 

(23e) 



As with the Mathieu eq., one can easily solve the unforced equations between kicks, where surface waves propagate 
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and are damped by viscosity. Let u n (z,t) be the z-component of the velocity valid for (nAi) < t < (n + (l/2)Ai). 
Let C n (f) be the amplitude of surface displacement over the same time period. From Eqs. (p3|), we find 

u n (z, t) = e^ (t-in) u n (.z) + c.c. u n (z) = a n e KZ + b n e kz (24a) 

C„(t) = Cne^-^ + c.c. k 2 = fc 2 + i/>/v (24b) 



Eqs. ( p4| ) give a n = —2vk 2 ( n a nd b n — (2vk 2 + ip)( n . 
The normal-stress condition (Eq. 23d ) gives the disper- 
sion relation ip(k) for the surface waves: 



ip = —2i>k 2 ± iojo\ 1 — 



4is 2 k 3 n 



(25) 



where ujq = gk + (j/p)k 3 . Note that this expression 
for ip resembles the one found for the Mathieu equation. 
(See Eq. ^ and the line following.) The laterally infi- 
nite container supports a continuum of modes indexed by 
the 2D wavenumber k = (k x , k y ). A mode of wavenum- 
ber k has a damping coefficient 7 w 2vk 2 . Damping 
gives a (second-order) decrease in the frequency, so that 
uj' {) = u)Q./l — 4,7/2 % 3k . One complication, relevant only for 



higher viscosities, is that Eq. ( |25| ) is an implicit equation 
for ip, since k = k(i/j). One can write (and solve) an 



explicit, fourth-order equation for i/j(k), but not much 
insight is gained. The important point is that ip{k) is 
the dispersion relation for free surface waves propagat- 
ing on a deep, viscous fluid. Such dispersion relations 
have been catalogued for a variety of physical situations 
(small or large depth, small or large viscosity, viscoelas- 
tic effects, surfactant layers on the surface, etc.). One 
can look up or calculate the appropriate tp(k) for each 
physical situation one wishes to study. The dispersion 
relation is completely independent of the way the waves 
are excited (parametric or direct forcing). For now, we 
leave il>{k) as arbitrary. 

The next step is to derive the matrices relating the 
complex amplitudes Cn+i to Cn+1/2 to £„. At (n+l/2)Ai, 
continuity of ( yields Re(C„ + i/ 2 ) = Re(C„e' / ' At / 2 ). The 
velocity-jump condition is derived by integrating the 
normal-stress condition Eq. (23d) over a small interval 
of time centered on (n + 1/2) At. We find 



Ax k 2 

(dzU n+1 / 2 \z=o)(t n+ i/ 2 ) - (d z u n \z=o)(t n +i/2) = H (- 1 ) 



C„e 0A */ 2 + c.c. 



where 



(dzUn+1/2) (tn+1/2) = Ka n+1/2 + kb n+1/2 + c.c. = [2^k 2 (k - k) + ^klCn+i/2 + c.c. 
Inserting this into Eq. (E6h, we find 



ImC n+1/2 = ImC„e^ At / 2 + ^ 



At (Im-0 — 2h>klmn) 



Re( n e 



ipAt/2 



(26) 



(27) 



(28) 



Define now ijj = —7 + iu)' and K = K r + i/tj. Then 
Cn+i = e^ At / 2 M 2 ( n+1/2 and C„+i/ 2 = e^^M^. 
Here, the matrices Mi and M 2 have exactly the same 
form as before, except that now a = 4:X pp k/[(At)(ui' — 
2vkKi)\. The tongues then are given by 2a(At) = 

cos i^pAt ±cosh^y a T^ difference is that in a laterally 

infinite container, all modes k are available. Thus, what- 
ever the driving frequency oj (or whatever the At), the 
response of all modes must be considered. This implies 
that the overall threshold e* discussed above is the rel- 
evant concept. Since the to = 2luq tongue is in general 
the easiest to excite, one expects waves at just half the 



driving frequency. 

For small damping (a = Avk 2 juj -c 1), one has 



7 ; 
of 



2vk 2 , k = k, 



2ujq, giving a threshold 



vku), 



(29) 



which is slightly higher than the corresponding condition 
for sine- wave forcing 0, a pp = Avk 2 uj. 
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IV. DRIVING BY AN ASYMMETRIC 
TRIANGLE WAVE 

Recently, Edwards and Fauve B have explored the ef- 
fects of exciting Faraday waves with a driving signal that 
is the sum of two sine waves of rationally related fre- 
quency. They showed that for fixed frequency ratio (5:4 
in their original work), the relative amplitude and phase 
of the two sine waves could be considered as independent 
control parameters. Varying these new parameters dra- 
matically affects the pattern-formation selection. For ex- 
ample, they produced hexagons and quasi-periodic pat- 
terns ( "quasi-patterns" ) under conditions where single- 
frequency forcing would lead to a stripe pattern. 

Here, we look at a simple model of two-frequency forc- 
ing, an asymmetric triangle- wave forcing (Fig. |l|c) . This 
waveform has the crucial feature that it is asymmetric 
under time reversal, f(t) ^ f(At — £), in contrast with 
the "single- frequency" driving waveform (Fig. [l]b). As 
Edwards and Fauve have argued, time-reversal invariance 
of f(t) implies that the nonlinear amplitude equations for 
the parametrically excited waves must be invariant un- 
der A — > —A; however, this symmetry is removed for 
asymmetric fit). 

In the delta- function model, asymmetric triangle-wave 
forcing is easily solved by replacing the At/2 in the ma- 
trices Mi and M 2 (Cf. Eqs. |o|) with At x and At 2 , where 
Aii + At 2 — At. The tongue boundaries for the Mathieu 
Eq. are then given by 



e c (At) = y/2u 



I cos uj' At ± cosh ^At 
sinc^QAii sinu^A^ 



(30) 



The tongues are plotted for several choices of At\ and 
At 2 in Fig. |c. Note the reappearance of the harmonic 
tongues for Ati ^ At 2 . In this simple model, there is 
no analog of the "bicritical point" discussed by Edwards 
and Fauve, where two different wavenumbers go unsta- 
ble simultaneously. Nonetheless, the symmetry argument 
given above suggests that asymmetric-triangle- wave forc- 
ing will alter the nonlinear pattern selection in much the 
same way as two frequencies do. 



V. CONCLUSION 

In this paper, we have shown how the linear-stability 
analysis of Faraday waves is greatly simplified for 
triangle-wave forcing. Our analysis is simple enough to 
include in an undergraduate hydrodynamics course. One 
could also combine the analysis with a simple experiment 
to make a good senior thesis project. References H and 
p5[ give careful discussions of the required experimental 
apparatus |Ts| ]. 

Our calculation may also serve as the starting point for 
a more extensive nonlinear analysis. To date, nonlinear 



amplitude equations for Faraday waves have been derived 
only in the limit of small damping a <C 1 ]l6|,[l7|]; how- 
ever, many interesting phenomena are seen when a ~ 0.4 
to 0.7 ||. With delta- function driving, the effects of the 
parametric pumping are straightforward, even in the non- 
linear regime. The hard part of the calculation becomes 
the determination of nonlinear terms in the dispersion 
relation for surface waves. However, this is a problem 
that has received much attention over the years, and a 
number of results are available [|19| . 
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